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This paper presents an integration of distributed generation and capacitor in 
radial distribution system via Firefly Algorithm (FA). In this study, the FA is 
developed in order to determine the optimal location and size for 
compensation schemes namely distributed generation (DG) and 
compensating capacitor (CC). The FA which is a meta-heuristic algorithm is 


inspired by the flashing behavior of fireflies. The proposed technique was 





tested on IEEE Reliability Test systems namely the IEEE 69-bus and the 
Keywords: program was developed using the MATLAB programming software. 
The results shown a significant reduction in the line losses and voltage 
profile improvement has been obtained with the installation of distributed 
generation and capacitor in the system. 
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1. INTRODUCTION 

For several years, a variety of techniques and procedures had been developed by many researchers 
for loss minimization issue and allocation of DG. This problem is a long progressing issue mainly in 
distribution system. However its effect led to temperature rise due to significant voltage drop through the 
distribution line. Thus, compensation scheme should be proper scheduled in the attempt to alleviate the 
voltage drop phenomenon. Voltage drop problem may cause insecure power delivery to consumer. 
The consequence of voltage drop is loss increment in the whole system. Thus it requires compensating 
scheme such as static Var compensation (SVC) [1], compensating capacitor [2] and DG [3-5]. In the present 
power system, DG installation has shown an escalating growth in the distribution networks around the world 
due to the raise in promotion towards utilization of renewable energy resources and development of co- 
generation plants. Distributed generations (DGs) is also known as embedded energy refers to power 
generation at the point of consumption. In other words, DGs is the process of generating electricity at a 
specific location and then connecting that supply into the electricity network. While capacitor bank is widely 
used in AC power supply for error correction. Generally, inductive load causes susceptible to phase shifts and 
power factor lags in the power supply cause a loss of system efficiency if left uncorrected. With proper 
planning, the integration of distributed generation and capacitor in a distribution system would lead to 
enhancement in the network performance in terms of the voltage profile improvement reduction in line losses 
and improve reliability and quality [6-9]. As the result, the total demand required from the grid could be 
reduced, thus cutting the need to strengthen the feeders connecting the network to the grid. 

Thus, the best location and sizing factors are importance issues in the implementation of distributed 
generation and capacitor bank in the distribution system[10]. With optimal placement and sizing of 
distribution and sizing of distributed generation or compensating capacitor or both of them, the maximum 
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possible advantages could be obtained especially to improve the performance and to reduce losses in the 
system [4, 11-14]. 

A lot of conventional analytical methods such as loss sensitivity method, Linear Programming (LP) 
[15-17] and Lagrangian approach [18, 19] have been applied to solve the allocation problem. Nevertheless, 
these methods are very time-consuming and low reliable, especially when dealing with complex power 
problems [20-22]. With the increase of computation intelligence technologies especially in Nature Inspired 
Algorithm (NIA), many researchers have turned their attention to meta-heuristic approach [23]. 
The integration of DG and CC units provides potential benefits to the conventional distribution systems. 

In this study, the new swarm intelligence approach that utilizes FA to determine the best location 
and sizing of DG, CC and both DG-CC in order to minimize the total system real power loss is proposed. 
By using FA, the size and location of proposed schemes can be determined simultaneously. This paper 
presents an integration of DG and CC in radial distribution system via Firefly Algorithm (FA) and the 
effectiveness of the proposed methodology was verified by the analysis on a 69-bus distribution system. 


2. RESEARCH METHOD 

The Firefly Algorithm (FA) is a metaheuristic nature-inspired optimization algorithm which it is 
invented by Yang for solving the multimodal optimization problem [24-26]. It creates based ib the behavior 
of fireflies, or lightning bugs. Fireflies have about more than two thousand species where the flashes often 
unique for each particular species. The flashing light is produced by a process of bioluminescence where it is 
functioning as signaling system. 

Generally, the flashes from firefly are to attract mating partners and to attract potential prey. Firefly 
algorithm has three particular ideal rules which are based on the behavior of firefly. These are following : 

a) All fireflies are unisexual so that any individual firefly will be attracted to all other fireflies; 

b) Attractiveness is proportional to their brightness, and for any two fireflies, the less bright one will be 
attracted by (and thus move towards) the brighter one; however, the intensity (apparent brightness) 
decrease as their mutual distance increases; 

c) If there are no fireflies brighter than a given firefly, it will move randomly. 


2.1. FA for Optimal Allocation and Sizing 

In this section, FA is proposed to find the optimal location and size of compensating schemes in the 
distribution system to minimize the total loss in the system. By minimizing the loss, the voltage profile at 
each bus is also expected to be improved. The variables of the optimal allocation of schemes are coded as 
follows: 


X= [yp X15) 00+ Xap Xns] (1) 


where L indicates the location of compensating scheme, S indicates the size of compensating scheme and n is 
the number of compensating scheme that need to be installed in the system. These variables then are included 
in the data of load flow and the load flow is executed to obtain the total loss of the system. 

The procedure to obtain the optimal allocation of schemes requires load flow to be run iteratively. 
After obtaining the best location and size of schemes simultaneously, the procedure is stopped. Results will 
be recorded. The objective function, f(x) is the results of the total loss of the system, Pros; to be minimized, 
as follows: 

f (x) = min(Xin Ploss) (2) 

Where line is number of transmission lines in the distribution system. 

The following procedures were implemented in the simulation in order to determine the optimal 
sizing of compensating scheme: 
Step 1: Generate control variable 
Initialization process was conducted by generating control variable using a uniformly distributed random 
number generation, X;. The random numbers represent the value of x; and 7; as the variable to be optimized. 
During initialization, the population size depends on the number of individuals. 
Step 2: Constraint violation test. 
Test the value of X;. The calculation of the fitness value (y;) of each individual is done to ensure the fitness is 
within the acceptable limit. 
Step 3: Selection process for ‘n’ fireflies 
Based on the ranking data, the best ‘n’ population is selected as initial location of ‘n’ numbers of fireflies. 
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Step 4: Calculate distance and attractiveness 

Move firefly i towards j in d-dimension. Attractiveness varies with distance r via exp '! Evaluate new 
solutions and update light intensity. The firefly attractiveness, 8 can be defined by (3). The distance between 
any two fireflies i and j are calculated using the Cartesian distance in Error! Reference source not found.. 
The movement of firefly i is attracted to another more attractive (brighter) firefly 7, determined by Error! 
Reference source not found.. 


- -" 














1 (4) 


2 
x, =x,+f ee (x .- x, )+a(rand - 0.5) 
i i 0 joi (5) 


Step 5: Fitness calculation 
Fitness calculation is conducted for all individuals in the population. Fitness value is the one which will 
undergo the optimization process. The fitness will be maximized or minimized depending on the objective 
function. 

The process of incorporating the FA into optimal allocation of compensating schemes is shown in 
Figure 1. 





Initialize location of fireflies 
Iteration=1 


Insert variable x into load flow data 
Run load flow 
Objective function evaluation 


Ranking fireflies by their light intensity/objective 


Find the current best solution 


Move all fireflies to the better locations (Updating 


fireflies) 


Iteration=50 


Yes 













Figure 1. Flowchart for FA algorithm 


3. RESULTS AND DISCUSSION 

The FA has been tested on IEEE 69-bus radial distribution system. Figure 2 shows the single line 
diagram for the test system with a total real and reactive power demand of 3802.19 kW and 2694.60 kVAr 
respectively. The FA properties in this simulation are set as follow: 


Number of fireflies: 20 Schemes size: 0.01 MW< Pychemes < 2.5 MW 
Maximum iteration: 30 Alpha (scaling parameter): 0.25 
Number of DG unit: 1 Minimum value of beta (attractiveness): 0.2 
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Number of CB unit: 1 Gamma (absorption coefficient): 1 
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Figure 2. IEEE 69-bus radial distribution system [27] 


The aim of this study is to see the impact of the proposed method on the IEEE 69 bus distribution 
system using different compensating sc. The following tables will show the result after carrying out the FA 
with random bus selection for 30 samples. Table 1 shows the results of objective function in (2) with the 
installation of single DG in the system. The minimum loss of 0.0232MW is experienced by the system when 
DG is installed at bus 61 with the loss reduction is recorded as 89.7% as highlighted in yellow. 


Table 1. Optimization with Single DG Placement 











No. Status DG - Total Losses Loss Reduction 
Location Sizing (MW) MW (%) 
1 65 1.3955 0.0617 72.56 
2 11 1.8082 0.1627 27.64 
3 15 0.8848 0.1823 18.95 
4 34 0.0364 0.2249 - 
5 24 0.7103 0.1882 16.33 
6 2 2.5 0.2249 0.01 
7 59 1.9118 0.0435 80.65 
8 2 2.5 0.2249 0.01 
9 15 0.8848 0.1823 18.95 
10 24 0.7103 0.1882 16.33 
11 20 0.7837 0.1851 17.71 
12 64 1.6076 0.0408 81.87 
13 20 0.7837 0.1851 17.71 
14 37 0.2146 0.2249 - 
15 64 1.6076 0.0408 81.87 
16 46 0.0976 0.2248 0.03 
17 53 2.5 0.1429 36.46 
18 41 0.1614 0.2249 0.01 
19 14 1.0038 0.18 19.98 
20 59 1.9118 0.0435 80.65 
21 40 0.1645 0.2249 0.00 
22 48 0.9295 0.2242 0.32 
23 53 2.5 0.1429 36.46 
25 31 0.0595 0.2249 - 
26 35 0.0275 0.2249 - 
27. 52 1.9695 0.174 22.63 
28 32 0.0526 0.2249 - 
29 40 0.1645 0.2249 0.00 
30 34 0.0364 0.2249 - 





Similar study was performed for the integration of single capacitor placement as compensating 
scheme. The minimum total loss occurs when DG is installed at bus 61 with the system loss of 0.152MW. 
The optimal capacitor size is 1.33MVar to achieve this result as tabulated in Table 2. 

Lastly, the determination of DG location for the installation of DG and capacitor simultaneously is 
also conducted. The results is tabulated in Table 3. The minimum loss of 0.0235 MW is experienced by the 
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system when DG is installed at bus 63 while capacitor at bus 24. This required a DG size of 0.0364 MW and 
a size of compensating capacitor is 0.3770 MVar to compensate the system for loss reduction. 
Table 2. Optimization with Single Capacitor Placement 

















No Status Capacitor Total Losses Loss Reduction 
. Location Sizing (MVAR) (MW) (%) 
1 62 1.3087 0.1529 32.00 
2 33 0.0307 0.2249 - 
3 27 0.3976 0.215 4.42 
4 18 0.5432 0.2119 5.77 
5 32 0.0332 0.2249 - 
6 10 1.3283 0.2045 9.06 
7 56 1.6687 0.1898 15.60 
8 57 1.4536 0.1729 23.12 
9 40 0.114 0.2249 - 
10 56 1.6687 0.1898 15.60 
12 30 0.0379 0.2249 - 
13 10 1.3284 0.2045 9.06 
14 17 0.5435 0.2119 5.78 
15 38 0.1182 0.2249 - 
16 62 1.3087 0.1529 32.00 
17 8 2.0722 0.2006 10.79 
18 41 0.0704 0.2249 - 
19 39 0.1143 0.2249 - 
20 60 1.359 0.157 30.18 
21 62 1.3087 0.1529 32.00 
22 64 1.1482 0.16 28.85 
23 37 0.1579 0.2249 - 
24 20 0.5094 0.2125 5.53 
25 12 1.008 0.2069 8.02 
26 55 1.7434 0.1927 14.30 
27 20 0.5094 0.2125 3,93 
28 6 2.1158 0.2137 4.97 
29 19 0.522 0.2123 5.62 
30 31 0.0366 0.2249 - 
Table 3. Optimization with DG and Capacitor Placement 
No Status DG Capacitor Total Losses Loss Reduction 
. Location Sizing (MW) __ Location Sizing (MVar) (MW) (%) 
21 1.3955 34 0.7609 0.1979 12.00 
54 1.8082 41 0.2161 0.1362 39.45 
66 0.8848 69 0.2660 0.1671 25.6 
40 0.7103 22 0.4650 0.2141 4.82 
52 25 66 0.7798 0.1679 25.33 
19 1.9118 54 1.4475 0.1649 26.70 
9 2D) 49 1.0083 0.1493 33.61 
5 0.8848 39 0.0534 0.2222 1.18 
2 0.7103 50 0.5279 0.2241 0.37 
17 0.7837 44 0.1611 0.1834 18.47 
52 1.6076 47 0.8475 0.1740 22.65 
39 0.7837 35 0.0014 0.2256 0.33 
4 0.2146 52 1.3795 0.2085 7.31 
3 1.6076 42 0.1230 0.2248 0.03 
65 0.0976 50 0.7580 0.0611 72.84 
40 2.5 8 1.9330 0.2007 10.75 
4 0.1614 42 0.0628 0.2247 0.11 
4 1.0038 19 0.5211 0.2120 5.15 
24 1.9118 6 1.5947 0.1814 19.34 
8 0.1645 21 0.2892 0.1498 33.41 
37 0.9295 56 1.6686 0.1898 15.61 
29 2,5 41 0.0900 0.2251 0.11 
51 1.8368 41 0.7308 0.1630 27.54 
62 0.0595 40 0.5212 0.0253 88.74 
51 0.0275 28 0.5070 0.1575 29.98 
58 1.9695 64 0.9180 0.0450 80.01 
21 0.0526 47 1.1209 0.1860 17.31 
44 0.1645 34 0.0239 0.2248 - 
67 0.0364 64 1.0300 0.1183 47.42 
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A comparative study was done in order to show the capability of the technique to perform DG 
planning and to ensure the proper placement and sizing of DG unit. The results show the percentage of 
reduction for installation of DG managed to outperform the other with considerable percentages of loss 
reduction of 89.70%. 


4. CONCLUSION 

This paper has presented the integration of DG and CC in Radial Distribution System via Firefly 
Algorithm has been successfully implemented and tested on JEEE 69-Bus Distribution System. 
The effectiveness of FA was demonstrated and tested. The results indicated that the proper planning of the 
integration of distributed generation and compensating capacitor in a distribution system would lead to 
enhancement in the network performance in terms of the loss reduction in the system. The result obtained 
from this study could be beneficial to the power system operator or to the utilities to identify the 
compensating scheme placement with a few selection locations to install the DG in the system. 
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